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In a recent paper we have analyzed the AdS / CFT duals to shock waves propagating in the 
M = 4 plasma. Here we study further properties of the system. In the gravity description we 
consider the properties of the dual black holes, showing in particular that they are stationary 
black holes with expanding horizons. This is possible because the horizon is not compact; 
in the fluid, this corresponds to the situation when entropy is being produced and carried 
away to infinity. We also consider shocks in dimensionalities d other than four and find that, 
for plasmas whose duals are given by asymptotically AdS spaces, the exponential tail of the 
shock on the supersonic side shrinks as 7~s as the velocity approaches the speed of light 
(the Lorentz factor 7 goes to infinity). This generalizes the behavior 7" 2 we have found 
previously for d = 4. Finally, we consider corrugations of the shock front and show that the 
shock is stable under such perturbations. There are, however, long lived modes, excitations 
of which describe generation of sound by the shock wave, the energy for this being provided 
by the incoming fluid. 
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1 Introduction and summary 



A useful tool for studying gauge theories at strong coupling is provided by the AdS/CFT 
correspondence [II El El E] • Lately much effort has been concentrated on studying the dy- 
namics of strongly coupled plasmas. Several reviews can be found in [SJ El El El El EE]- In 
conformal plasmas at strong coupling the hydrodynamic description is valid up to length 
scales of order of the inverse temperature. At shorter scales, hydrodynamics is not valid and 
one must study the microscopies of the theory. For theories with AdS duals, a way beyond 
the hydrodynamic description is to study the dual gravity theory, an approximation which 
does not break until the much smaller string length scale. Within that context we wish to 
study a specific phenomenon, shock waves in conformal plasmas. 

Shock waves play an important role in the dynamic of fluids in the supersonic regime. 
In recent years, a particular fluid, the quark gluon plasma, has been studied experimentally 
and evidence has been found for creation of a Mach cone by a moving heavy quark. This 
QCD plasma is difficult to study theoretically because it is in the regime of strong coupling. 
One intriguing possibility to gain some insight into its properties is to use the AdS/CFT 
correspondence, which maps the physics of a closely related plasma to the physics of black 
holes in AdS space. The plasma in question is a conformal plasma that appears when the 
M = 4 SYM is heated to a finite temperature. Although not the same as the QCD plasma, 
it is strongly coupled and may have similar properties, at least at the qualitative level. 

Shock waves can appear during the initial stages of creation and thermalization of the 
Quark Gluon plasma. They create entropy and might be responsible for the whole initial 
entropy creation; in certain models, as for example in [llj, the latter determines the number 
of produced particles. Shocks can also appear when a hard parton moves with a supersonic 
velocity. In this case, one expects that they are closely related to the parton's Mach cone. 
Several aspects of this process have been studied within AdS/CFT [121 131 EH EEl EE]- In 
the heavy ion collisions literature the existence of shock waves and their consequences have 
been studied both theoretically and experimentally [IZl HH HH [20l EH [221 EH EH [251 ESI [271 

ESI EH E01 EH E21 E31 EU- 

In a recent paper [35], we have observed that, since shock waves probe the microscopic of 
the theory, they are in the regime where, according to AdS/CFT, the dual gravity description 
becomes important. More precisely, for a perturbation of wavelength much larger that 1/T, 
where T is the temperature, we can just use hydrodynamics but for shorter wavelengths the 
appropriate microscopic description is the gravity dual. In fact, the correspondence allows us 
to work both ways: when considering weak shocks, which can be resolved in hydrodynamics, 



2 



we have learned about the existence of dual black holes with solitary waves propagating on 
their horizons, while when considering strong shocks we have been able to use gravity to 
compute the penetration depth of the exponential tails on both sides of the shock and thus 
learn something about the fluid. 

In this paper we continue our investigation of these issues. We study in more detail 
the gravity solutions corresponding to weak shocks, especially, in regard to the properties 
of their horizon and the fact that they create entropy even when they are stationary. For 
strong shocks, we study the properties of the exponential tails in dimensionalities other than 
four and show that, just as in d = 4, for shocks propagating close to the speed of light 
(7 — > 00), the scaling of the penetration depth on the supersonic side is determined by the 
near-boundary region of the dual geometry and, in this sense, is universal for all fluids with 
asymptotically AdS duals. We also explore a region close to the horizon where there are 
potential divergences. Finally, we investigate the question of corrugation stability of the 
shock and show that the shock can produce sound, which travels in a definite direction that 
depends only on the velocity of the shock. 

2 Shock waves in ideal and first order hydrodynamics 
2.1 Shock waves in ideal hydrodynamics 

In this section we examine shock waves in ideal hydrodynamics in arbitrary dimensions. The 
stress energy tensor for an ideal conformal fluid in (d — 1) space and one time dimension can 
be written as 



where p is the pressure, T the temperature, the <i-velocity field, and G N is a constant 
that in holographic computations is naturally identified with Newton's constant in the dual 
d + 1 dimensional gravity model. We follow the notation of [36J which arises from the 
holographic computation of ([!]). For reference we give the value of in three interesting 

cases. The most well studied AdS/CFT duality is between d — 4, J\f — 4 SU(N) super 
Yang-Mills theory and type I IB string theory on AdS§ x S 5 . In this case G^ = . For 
d — 3, J\f — 6 super conformal U(N)k x U(N)_k Chern-Simons theory is believed to be dual 
to M-theory on AdS 5 x S 5 [37], with G$ = ^ where A = f is the 't Hooft coupling. 
For the d = 6 superconformal theory living on N M5 branes we have that G^j = 5^3- It 
should be noted that in the following whenever we do not explicitly give the dependence on 




p = 
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x = 



Figure 1: Sketch of a shock wave in the rest frame of the interface. For ideal hydrodynamics 
there is a discontinuity at x — 0. Including the higher order terms in the expression of the 
stress energy tensor resolves the discontinuity. The conventions of the paper are that the 
fluid moves to the right. The left hand side of the flow is supersonic and the right hand side 
is subsonic, the choice prescribed by the second law of thermodynamics. 

G^ +1 ^ we have normalized the stress energy tensor by 167tG^ +1 \ 

As in [35] we study a shock moving in the x direction in the frame in which the shock 
front is static. The two components of the stress tensor that are constant are T tx and T xx . 
In the present case, these are given by 



1 - r- " " 1 -t; 2 



T^=pd T ^-=, 2W =p (i + d -!!—). (2) 



Labelling the fields on the supersonic side with subscript 1 and those on the subsonic side 
with subscript 2, we find that the fields are related by 



2„,2 1 rp / ( A -i \2„,2 1 \ l / d 

(3) 



1 p 2 {d-lfvl-l T 2 f{d-lfvl-l 

V\V2 



vi ' T x \{d - 1)(1 - v\ 

For d — 4 these agree with the expressions derived in [35]. It should be noted here that 
we study shocks for d > 2. The case d — 2 is special since the speed of sound coincides 
with the speed of light. There are no supersonic flows for conformal fluids in two spacetime 
dimensions. This is slightly disappointing: as noted in [SB] conformal fluid dynamics in two 
dimensions is trivial Q There are no non-zero Weyl covariant tensors that can be constructed 
from the velocity and temperature fields. 

The shocks are the only mechanism that generates entropy in ideal hydrodynamics. The 



1 Since left and right moving waves do not interact in two dimensions there is no local equilibration process. 
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difference in the entropy ffux between the two regions is given by 

47rT 1 \ d ~ 1 / Vl {{d-lfvl-lf'' 2 



d 



2d 



8^-8^ = 4* — ^=^"^ 71 ~h=T ■ ( 4 ) 



((d-l)(l -«?))- 



2.2 Shock waves in first order hydrodynamics 

When we include effects of viscosity, the profile of the shock is no longer discontinuous. The 
temperature and velocity fields are varying continuously. The first correction to the ideal 
stress energy tensor ([T| is given by [Ml E] 

= ^ + u,u u , a*"' = P» a P^d [a u p) - J^P^d a u a . (6) 

The equations of motion are given by d^T^ u = 0, where the stress energy tensor is now the 
sum of (jl]) and (|5]). The nontrivial equations are those for the T tx and T xx components. To 
obtain these components, we first compute 



d-2 . . , r -^z , , , d-2 , 



-u x (x)y/l + u x (x) 2 u' x (x) , a xx = -r—^ x {x) (l + u x (x) 2 ) . (7) 



d 

Then, the components of the stress energy tensor are 

d-l 



AvrTV 1 . . , — / d-2 ,. A 

T tx = I —t- I u x [x)y/l + u x {x) \AixT - 2-—jU x (x) I 



47rT\ d ~ 1 /4ttT, . , d 



1 + du x (x) 2 ) - 2— — (1 + u x (x) 2 )u' x (x) ) . (9) 



d J \ d d — 1 

We expect the first order hydrodynamics to be a good approximation for perturbations 
whose characteristic wave numbers in the rest frame of the fluid are small compared to the 
temperature: 

|g'|= 7 |g|<T. (10) 

It turns out that for shock waves this condition is satisfied only when the velocities on both 
sides of the shock are close to the speed of sound. We then search for a solution for the x 
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component of the rf-velocity field in the form 



u x (x) = y -^-^ + u 00 Su(x) , (11) 

where Uoo = u x (—oo) is a small parameter. It is convenient to introduce, instead of x, a 
rescaled variable 

* = d=l X ■ (12) 

To the leading nontrivial order in u^, the equation of motion for 5u(£) becomes 

Su'(0 = 5u(0 2 - 1 , (13) 

which has the following solution: 

5u(0 = tanh(-0 . (14) 

The temperature perturbation can then be obtained by noting that the stress energy tensor 
satisfies the "Landau gauge" condition 



(4irT) d 

or more explicitly 



T^ = -{d-l)^^u U) (15) 



7 / \ l/d 

T= , , , ( T xx_ T tx^j\ _ ( 16 ) 



4n(d-l) l / d V u x J 
Asymptotically, the behavior of the solution is exponentially decaying 



, d-2 

with 



" i (•'• ) - \ / -t— ^ TUoo + A±e i *T x , (17) 



lS_ = J d _^l v = A —u (18) 



where 



d-l^ 3/2 



(19) 



d-2 / 

Another way to use the first-order hydrodynamics is to perform linearized analysis around 
a fluid moving uniformly at constant temperature. Formally, this can be done for any fluid 
velocity, but as we will see in a moment it is meaningful only for velocities close to the speed 
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of sound. Assuming the following form for the velocity and temperature fields 



u(x) 



v 



— V 2 



iqx 

+ due^r 



x) 



iqx 

T + 6Te^ 



(20) 



and linearizing the equations of motion, one obtains a 2 x 2 linear system for Su and ST. For 
this system to have a nontrivial solution the determinant must be zero and this determines 
q to be 

iq 

7rT ~ (d - 2)v 



2VT 



{v 2 (d-l)-l) 



(21) 



For v close to the speed of sound, \q\ is small, and (21) agrees with the result (18) obtained 



from the solution to the nonlinear problem. For other values of v, the small-gradient condi- 



tion (10) breaks down and the first-order theory is inapplicable. Note in particular that for 



ultrarrelativistic speeds (v — >■ 1) q as given by (21) tends to zero in disagreement with the 



result that we obtain from the analysis of the black hole quasinormal modes in Section [3j 

The entropy current in the first order hydrodynamics is given by s M — 47r ( 47rT ) 

and satisfies 



16ttG 



id+i] 



2r] 



fW 



2r] d-2 



u' x (xf 



(22) 



T ^ T(x)d-1 

Thus, the entropy production in the first order theory is symmetric around the position of 
the shock front (x - 



3 Linear analysis for strong shocks 

Strong shocks are characterized by large gradients of the velocity and temperature. As such, 
they cannot be analyzed within hydrodynamics. Since the gauge/gravity correspondence 
is not limited to large wavelengths, it can provide information about strong shocks. In 
particular, we can carry out linearized analysis near a uniform flow directly in the dual 
gravitational theory. 

Consider then a uniformly boosted black hole in d + 1 dimensions. The metric is 

, 2 dz 2 z d ~ 2 , , n , . , „ , x 2 -dt 2 + dx 2 

ds = ~T( —v + —r(™shpdt-smhpdx) + . (23) 

z 2 (l-(^) d ) Z h z 

For convenience we have set L, the radius of AdS equal to 1. The parameter /3 is related to 

2 In the next (second) order hydrodynamics, there is an ambiguity in the definition and divergence of the 
entropy current. For a relevant discussion one can see [39] . 



7 



the velocity of the boost v via 

v = tanh/3. (24) 

The temperature of the dual plasma can be computed from the location of the horizon Zh 
and we find that 

T=j d ~. (25) 



Via the standard holographic dictionary [21 E], the metric (23) gives the stress energy tensor 
([I]) of an ideal fluid moving with velocity v. The normalization of ([!]) has been chosen so 
that it naturally arises from the gravity computation. We are interested in perturbations in 
the sound channel, which are given by 

e iqx 

ds lert = — a" { H oodt 2 + H u dx 2 + 2H m dtdx + Hdx 2 d _ 2 ) . (26) 

Interestingly, in the case of the sound mode the following linear combination decouples and 
can be studied separately: 

Z(z) = H m {z) + g{z)H{z) , g(z) = 1 + ^j^~ 2 ■ (27) 

i 1 — v z 

We are considering here a perturbation with no time dependence: u> = 0, that is we are 
asking what is the penetration length for a static perturbation around a boosted black hole. 
If we think of that perturbation as being due to a shock wave present somewhere in the fluid, 
the penetration length can be interpreted as the width of the exponential tail of the shock. 

The linearized equation of motion for Z(z) can be derived from Einstein's equations as 
explained in appendix (TaI) and is 



Z"{z) + P{z)Z'{z) + Q{z)Z{z) = , (28) 

where 

-{2d - l)z" + z d (3d-l- 2^(1 - - 2jgg(l - tg) 
[Z) z(^-l)(^-2H(l-^) ' 1 1 

d 2 z 2d ~ 2 2 z d + 1 - v 2 

Q[Z) (z d - 1) {z d - 2^(1 - v*)) + q (z d - 1) 2 (1 - v*) ■ (30) 



In (27)-(30) we have set zy t = 1. The full dependence on the location of the horizon z^ can 



be recovered by the rescaling 

z 

z -> — , q-+ qz h . (31) 

Zh 



S 



3.0.1 Behavior close to the pole 

The equations of motion for the sound mode have a pole at 



d-2 

The pole coincides with the location of the horizon when 



z ' = 1 ^"' = ^W^T)- (33) 

This velocity is always greater than the speed of sound. For velocities greater than Vf the 
pole lies between the boundary and the horizon, the physical region for which we wish to 



solve (28). Therefore it is important to understand the behavior of Z(z) close to the pole. 



Expanding (28) around Zf we find that there are two solutions 

Z(z) ~ (z - z f ) s , s = , s = 3. (34) 

Fuch's theorem [40] guarantees that the larger root s = 3 corresponds to a regular solution. 
The second solution can in general have the form: 

Z(z) ~J2 Cn ( z ~ + H~ Cn + ? ^ z - z fY + ^°?>( z - Zf) , (35) 

n=0 n=0 

that is there can be logarithmic terms starting with {z — Zf) 3 log(z — Zf). Whether such log- 
arithms are in fact present depends on the precise values of the coeffcients in the expansions 
of P{z) and Q(z): 

oo 

p ( z ) = ^ZY + ^ Pn{z ' Zf)n ' (36) 

Z Z f n=0 

oo 

w = jzj +E^( 2 : - z f) n ■ ( 37 ) 

Z Z f n =0 



Substituting (35) into the equation (28), we obtain a recursion relation for the coefficients 



c n , c n . Specifically, cancellation of the (z — zf) 1 poles gives 



zt 1 



f 

c l = c 0^7 7T > ( 38 ) 



2(z d f -r 
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cancellation of the (z — Zf)° terms gives 



C2 = C 



2q 2 + z d f - 2 (d- l)(d-2) 
A(z d -l)(d-2) 



d 



(39) 



while cancellation of the (z — Zf) 1 terms determines 63. For C3 to vanish, the coefficients of 



(36), (37) have to satisfy 



P-i = -2 , (40) 
2 (p <?o + qi) + qi {vl + Pi) = • (41) 

These relations are indeed satisfied for all values of d. As a result, no logarithmic terms 
appear, and both solutions are regular at z — Zf. 

It should also be noted that the coefficients C\,C2 diverge as Zf approaches 1. The case 



Zf = 1 has therefore to be examined separately. When (33) is satisfied (28) becomes 



Z"\z] + 



d-l + z d {l - 2d) 



Z'(z) + 



2^2(2-2 



d 2 Z 



,2- d + 2z d (d-l) 



z{z d -l) v ' \(z d -l) 2 * (d-2)(z d -l) 2 

The solution close to the pole is Z(z) ~ (1 — z) s with 



Z(z) = . (42) 



s = 1± 



iq 



y/d(d - 2) 



(43) 



It is then convenient to define a function $ that has a regular Taylor expansion near z — 1, 
via 



Z(z) = (1 - z) s <S>(z). 



Next, one observes that for 



d(d-2) 



(44) 

(45) 
4 

the coefficients in the equation for $ are all regular at z — 1. Therefore, both iq = 4- V d ^f 2 - 
are eigenvalues. However, since Vf is supersonic, the perturbation must decrease toward 
negative x and so only ^q = is accepted as physical. 
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3.0.2 Boundary conditions 

At z — y 0, the spacetime becomes asymptotically AdS with the boundary at z — 0. Close to 
the boundary the two independent solutions behave as 

Z x {z) rsj constant , Z 2 (z) ~ z d . (46) 

The first solution will change the metric on the boundary Since we do not wish that, we 
impose the boundary condition 

Z(0) = . (47) 

Following the holographic renormalization procedure jH] we also find that the boundary 
stress energy tensor has components 



(46 



Close to the horizon the two independent solutions are the ingoing and outgoing waves: 

Z in (*)~(l-^v&, Z^)~(l-,pA. (49) 

The terminology comes about because of the way the corresponding (time-dependent) per- 
turbations behave in the unboosted black brane background. Note that we need these 
expressions at imaginary values of q (as appropriate for perturbations decaying in space), 
so an analytical continuation is implied. We impose the infalling boundary condition at 



the horizon (that is choose the first of the asymptotics (49)). For the specifics of whether 
the chosen boundary conditions set a well posed problem, an interested reader can consult 
section (4.2) of [35J. 

We set up the numerical problem by looking for the solution as a Taylor expansion near 



the boundary and, after peeling off the singular part as in (44), a Taylor expansion near the 
horizon. Matching the logarithmic derivatives of these two expansions at an intermediate 
point determines the eigenvalue q. The coefficients of the expansion near the horizon grow 



when the spurious pole Zf approaches 1 (cf. (38) and (39))), which requires us to keep a 
large number of terms in the expansion (we go up to around thirty). This procedure can 
be further checked in two ways. First, since we have seen that the solution is regular at 
Zf, we can develop a third expansion near z — Zf. We then match this expansion with the 
expansions near z = and z = 1 at two intermediate points and find that the result coincides 
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to five significant digits with result from the previous procedure. Second, we can use the 



shooting method (integrate numerically (28) from the horizon to the boundary). For this, 



we circumvent the spurious pole by displacing it into the complex plane (zf —*Zf + ie). For 
a sufficiently small e, the results agree well with those of the other two methods. The three 
methods we have used are summarized graphically in fig ures [2(iQp(b)p(c]| 



3.1 Quasinormal modes for special values of the boost velocity 
3.1.1 Speed of sound 



When v 
are 



J-r, (28) simplifies and q = is an eigenvalue. The two independent solutions 

(50) 



Z(z) = z d , Z(z) =4 + z d log(l - z d ) 
Only the first one satisfies the boundary conditions. 



3.1.2 v 



For v = Vf 



2(d-l) 



j—^ the locations of the horizon and of the spurious pole coincide. The 



2(d-l) 



eigenvalue for this case has been found in sec. 3.0.1 



iq 



^d(d - 2) 



(51) 



It is interesting to note that, for this value of q, the boundary exponent (43 ), valid for the case 



zj = 1, coincides with the exponent in (49), derived originally for zj ^ 1. One consequence 
of this is that the curve q(v) is continuous at v = VfQ 

For d = 3, the wavefunction Z(z) corresponding to this value of q can be written in terms 
of a hypergeometric function: 



Z(z) = C 1 Z 3dl + C 2 Z : 



3(22 



where 



Z 3d i(z) = (z + 2)(z 2 + z + 1) 1/4 Vz^l 



2iV3z + iVS + 3 
3 — i\/3 — 2i\/3z 



iV3/4 



(52) 



(53) 



3 We are talking here about the "main" branch of q(v), which corresponds to the slowest-decaying com- 
ponent of a perturbation. This is the branch we are primarily interested in, although occasionally we will 
also present results for faster-decaying branches. 
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Im(z) 




(a) The Taylor series around the bound- 
ary z = and the horizon z — 1 are 
matched at an intermediate point. The 
location of the spurious pole z — Zf is not 
taken into account. When Zf > 1 there is 
no problem as the spurious pole does not 
lie in the physical region. When Zf < 1 
the spurious pole is in the physical region. 
Interestingly it does not affect the result 
of the numerical matching as long as it is 
not very close to the horizon. When the 
two poles coincide a new Taylor series has 
to be developed. 

A Im(z) 



A Re(z) 



z = 1 




(b) When the spurious pole is in the 
physical region we can develop three dif- 
ferent series, around z = 0, z = Zf, z = 1. 
The series around z — Zf has two free pa- 
rameters since there is no boundary con- 
dition to take into account. It is matched 
with the other two series at intermedi- 
ate points z\ (0 < z\ < Zf) and Z2 
[zt < z-i < 1). The result for q is numer- 
ically indistinguishable from the method 
in p(a) |. 



2 = 2/ + i( 




2 = 1 



(c) A third method used for determining 
q involves expanding around z = 1 to a 
low order in Taylor series and then in- 
tegrating numerically from the horizon. 
The spurious pole has to be displaced in 
the complex plane for the numerical inte- 
gration to be reliable. The sign of e has 
no effect and, for sufficiently small |e|, the 
results agree well with those of the previ- 
ous methods. 

Figure 2: Different methods for determining q 
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6 



4(>/3 + i) \/ C* 2 + 2 + l)(2iV3s + n/3 + 3) 
y/% — iV% — 2iV3z 

3(z + 2) (l + n/3 - 2z(l -iV3 + z)) l-iVg 3-zVg 6? 

+ z* + z + l 2i<l(1) 2 ' ; ~ 3i + v^3 + 2V3 2 J 

_ (54) 

3e 2^3 

^2 = ^1 7 7 777^- (55) 

Z — d 2-TlU, 2 ' 2 ' 2 ) 

It can be directly checked that this wavefunction satisfies both boundary conditions. 
In four dimensions Z(z) was determined in [35] to be 



Z(z) = -Re(y(z)) + C 1 Im{y{z)) , (56) 

= *VT^i (— j 2F1 (^-, l + ^■^ + - 2 y (57) 

where the constant C\ is numerically determined to be C\ = 0.38898. In six dimensions we 
were not able to find a simple wavefunction. 

3.1.3 The ultrarelativistic regime 7—7-00 

When the boost velocity approaches the speed of light the location of the spurious pole 
approaches the boundary. An estimate for q(v) in this case can obtained by using a small 



number of terms in the (35) and imposing the boundary condition Z(0) = 0. Retaining only 



the first three terms (corresponding to the coefficients Co, ci, C2), we obtain the estimate 

«L = L ^ ■ (58) 



Substituting the value of Zf from (32) and retaining only the leading terms in the limit 
v — > 1, we obtain 



The higher-order terms in the expansion (35) are non-negligible and significantly change 
the numerical coefficient in front of r ) 2 / d . In what follows, we determine that coefficient 
numerically. However, the 7 2//d scaling given by this simple argument is correct. 
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A more rigorous way to obtain this scaling is the following. Numerical solution of (28) 
suggests that when v — > 1 the eigenvalues q become large and the maxima of the wavefunc- 
tions are close to z = zj — > 0. One then is compelled to ignore factors of z compared to 



unity in the equation of motion for Z(z) (28). This changes the functions P(z) } Q(z) to 



_ 2(d-l) 2 - 7 V(3^-7d + 2) 
[) z{z d j%d-2)-2(d-l)) ' 1 ' 

Oto = - wf^'^oC 2 ,^ + - 1)) • (61) 



{z d 1 2 {d - 2) - 2{d - 1)) 
We can now define a new variable 



z 2/d z{d-2f/ d , 
x = — = 7 2 / d ^ ^ , 62 



and take the limit 7 — >• 00 keeping x fixed. Equation (28) transforms into 



Z"(x) - Z'(x) {3d d + Z(x)p 2 (-1 + 2 d —\x d ) = . (63) 



x[x d - 1) \ d - 2 

where 



p 2 = zjq 2 . (64) 



Nothing in (63) depends on 7. The eigenvalue p for this problem is therefore a simple 



number. The scaling for q then immediately follows from (64) 



1 2 

g ~ — ~1<i f (7 -» cx)) . (65) 
z f 

The physical region < z < 1 gets mapped to < x < 00 with Dirichlet conditions at both 
ends. The spurious pole at x = 1 can be circumvented in numerical integration by moving 
it in the complex plane to x = 1 + ie and then taking the limit e — > 0. Numerically, we find 
the following for physically interesting dimensions: 

iq = 1.019 7 2/3 , d = 3 , (66) 

ig = 1.8957 1/2 , d = 4 , (67) 
zg = 3.62 7 1/3 , d = 6 . (68) 
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3.2 Quasinormal modes for the boosted black hole in AdS± 

One physically interesting case is d — 3. The relevant SCFT is the Chern-Simons theory 
described by [37] and referred in the literature as ABJM, heated to a temperature T. The 
transport coefficients up to the second order have been calculated in (32] ■ Here, we find the 
values of q(v) for various v in the region < v < 1. The highest 7 factor for which we 
determine q is 7 = 35. Some interesting properties of the q(v) curve are as follows. 

For t! = 0we find iq = —2.1486. This is the mass of the lowest-mass glueball in QCD 2 Q 
It should be pointed out, though, that there are no shock wave solutions for which v = 
(on either side), so this quasinormal mode is not directly relevant to our problem. 

For v = 1/a/2 (the speed of sound) we find iq = 0, as expected. For small departures 
from the speed of sound, v ~ l/y/2 + 5v , we find, to the linear order in 5v, 

iq(v) = 4.24264 5v = 3^/2 5v . (69) 



This is in agreement with the hydrodynamic result of section (2.2) (provided we rescale (69) 
by a factor = ^j- ). For v = v j = we find numerically iq = 0.86602 = y/3/2 in perfect 



agreement with our analytical result (51). We also numerically explore the relativistic limit 
and find agreement with the analytical computations of Section |3. 1 .3 The numerical curve 
q(v) is shown in figure (j3]). 

It should be noted here that a simple approximation reproduces the curve quantitatively 
up to a 5% accuracy, with the largest error coming from the hydrodynamic region where iq 
is small. The approximation is given by 

2/3 Av + B , . 

lq = ^F^G' (70) 

where 

A=l, B = -l/V2, F = 0.1539, (3 = 1/(2.1486^). (71) 



The values of A,B,F,G are determined by requiring that the curve matches the numerics 
at the three points v — 0, 4g, w|. The function (70) does not change under the rescaling of 



all the coefficients so only three constants are independent. If one boosts back to the frame 
where the shock front is moving, this curve can give a good approximation for the location 
of the quasinormal modes on the complex u plane when u/ = — vq' . Such a determination 



4 We define QCD2 as a theory obtained by a Scherk Schwarz compactification [43], [44] of ABJM down to 
1 + 1 dimensions. 
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might be relevant if an effective theory of hydrodynamics of the type described in (13] needs 
to be pursued. Other values of A, B, F, G can be obtained if the point v = is replaced 
by the asymptotic behavior at v — >■ 1. This produces more accurate results for the region 

v>v f = \f\- 

3.3 Quasinormal modes for the boosted black hole in AdS? 

Another physically interesting case is d = 6. This theory is dual to a stack of N M5 
branes. The transport coefficients for this theory were obtained in [461 136] as a part of the 
calculation for generic dimension. The numerical method we have used to determine q is 
the expansion in Taylor series around the boundary and the horizon and matching of the 
logarithmic derivatives at a point between and 1. The numerical result up to 7 = 10 is 



shown in figure (3.3).. 

For v = 0, we find that the lowest glueball mass for the dimensionally reduced theory is 
given by iq = —2.703. For the speed of sound (v = l/\/5) iq = 0, and for small deviations 
v = 1/ yfb + Sv we find 

iq = 6.7082 Sv = 3^ Sv , (72) 



in agreement with (21 ) (remembering that Zh = in this case). For the case of the spurious 



pole coinciding with the horizon the numerical result agrees with the analytically predicted 



one. Finally the relativistic regime is found to be well described by (51). 



There is a simple approximation for the q(v) curve similar to (70): 



1/3 Av + B 

iq = 1 FvTG' (73) 

where now 

A = 1 , B = -l/VE , F = -0.0125 , G = 1/(2.7034^) . (74) 



4 Properties of the second-order metric 
4.1 Expression for the metric 

We now turn to the case of weak shocks, where we would like to use gradient expansion to 
infer properties of black holes in the dual gravity theory. In this section, as well as the next 
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(a) Values of Im(q) versus the velocity. The blue curve is the main 
branch, the green the hrst upper branch, and the dashed red the 
hydrodynamical approximation. In the hydrodynamical approxi- 
mation Im(q) diverges at v = 0. It is only a good approximation 

close to v — v s = v/i 



-Im(q) 




y 



(b) Values of Im(q) versus 7 = —^=5- Both the main (blue) and 

the upper (green) branch diverge at large 7 as 7 2 / 3 . We have de- 
termined Im(q) beyond the values shown here up to 7 = 35. The 
hydrodynamical approximation (red) gives vanishing q for large 7. 

Figure 3: Numerical values of q for different velocities in d = 3. 
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(a) Values of Im(q) versus the velocity for d — 6. The blue curve is 
the main branch and the dashed red the hydrodynamical approx- 
imation. In the hydrodynamical approximation Im(q) diverges at 

v = 0. It is only a good approximation close to v = v s = *J\- We 
do not show here the upper branch. 



-Ira(q) 




-2 



(b) Values of Im(q) versus 7 = ~7j=f ■ The main (blue) branch 

diverges as 7 1 / 3 . The hydrodynamical approximation (red) gives 
vanishing q for large 7 

Figure 4: Numerical values of q for different velocities in d — 6 
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one, we restrict ourselves to the case d = 4. 

For any solution of fluid mechanics that has small gradients of the fluid velocity and 
temperature, the dual bulk metric can be reconstructed order by order in the gradient 
expansion [47] . In our previous paper [35], we have found the first few terms in this expansion 
for planar weak shocks in the M = 4 fluid that extends indefinitely in the direction of 
propagation (the x direction). The expansion parameter is in this case Moo, the first-order 



variation of the asymptotic 4-velocities relative to the speed of sound (cf. sec. (2.2)). This 
parameter controls the size of the gradient terms. 

Here, we discuss some properties of the metric computed to the second order in M m . It 
has the following structure: 



9mn[t,x) 



( g tt (r,x) u°(x) g tx (r,x) 

u°(x) — u(x) 

9tx(r,x) -u(x) g xx (t,x) 

\ 













r 2 — \st{t, x) 

r 2 — x) j 



(75) 



We denote by u the component of 4-velocity of the fluid (relative to the shock) in the x 



direction and by u° the corresponding temporal component u° = \/l + u 2 . The function 
ST{r, x), which appears in the transverse components of the metric, is of the second order in 
and is given by eq. (84) below. 



Each of the metric functions in (75) includes a piece that has the same form as the 
equilibrium metric g^ N but now with slowly varying parameters, plus a "standalone" second- 
order correction: 

9MN(r,x) = g { ^ N [r,u(x),T(x)\ + H MN (r,x) . 
Explicitly, the equilibrium metric is 



/ 



2 , 4(« ) 2 



u 

4 



U 





-u 













— u 

4 












\ 





r 2 J 



(76) 
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with 7Y = vrT. The slowly-varying parameters (to the second order) are 

u(x) = +u {1) (x) +u {2) (x) , (77) 
T(x) = T {0) ll-^u m (x))+T {2) (x), (78) 



47rT( )M oo X 



where 

U(i)(x) = -Moo tanh ^ , (79) 

while U{2) and T( 2 \ are related as follows: 

nT {2) = - ^7rT (0) u (2) . (80) 

Explicit expressions for «( 2 ) an d 7} 2 ) individually are also available [35], but will not be 
needed in what follows. 

Note that um depends on two constants, Tfo) and Mqo, the first of which can be thought 
as the "base" temperature of the fluid, and the second as the strength of the shock. The 
gradient expansion, in the present case, is an expansion in powers of Uoo <C 1. The value 



l/y/2 for the leading term in (77) reflects the fact that weak shocks are only possible for 
velocities close to the speed of sound. 

Finally, the "standalone" corrections are 

H tt (r,x) = -s T (r,x) , (81) 

1 \/3 

H tx (r,x) = —ru' (1) (x)-—s T (r,x) } (82) 

2 3 

H xx (r,x) = --ru {1) (x) + -s T {r,x) . (83) 



where 



with 



*r(r,x) = VZ^ F i u '(i)( x ) (84) 

(l + y2)(l + y)2 

Fi(y) =ln a -2arctan?/ + 7r, (85) 

y 



and primes denote derivatives with respect to x. 
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4.2 Location of the horizon 

Define the horizon of the second-order metric ( 75 ) as the surface 



r = th(x) 



(86) 



at which the normal vector 



^ = (0,1,-7^,0,0) 



becomes null. This is an apparent horizon defined similarly to the ones in 
expect that, just as the metric itself, r# is a sum of the slowly varying 

tt{x) = 7iT(x) 



(87) 
. We 

(88) 



and a standalone second-order term. The standalone term is sourced by r' T and so is 0(u^). 
In addition, it varies slowly with x, so its derivative is O(w^). Thus, to the second order 
order inclusive 

r' H ( x ) = r' T (x) = — — 7rT (0 )u' (1) (a;) . 



The null condition becomes 



g rr (x) - 2ur' T (x) = 



(89) 
(90) 



where to the required accuracy we can set u — 1/ a/2. 
A bit of algebra yields 



g (r, x) = r 



- 3 ra (i)W ; 



(91) 



for r — > rn(x) (and to the second order in Woo)- When this is used in (90), there is a curious 



cancellation between the last terms in (90) and (91). As a result, the equality r# ~ is 
good to the second order inclusive, i.e., 



r H (x) = 7iT(x) + 0(ul) . 



(92) 



In other words, at the second order, r#(x) tracks point by point the slowly varying temper- 
ature (78). 
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4.3 Killing vectors 



The second-order metric (75) has Killing vectors corresponding to translations in t, y, and 



z. The zeroth-order metric, given by (76) with ^-independent parameters, has in addition 
a Killing vector corresponding to translations in x. For that zero-order metric, a linear 
combination of d x and d t with appropriate coefficients yields the Killing vector 

X fiM = ( M °,0, M ,0,0), (93) 



for which the horizon of (76) at r = is a Killing horizon, i.e., xfo) is normal to it at every 
point. 

In general, existence of a Killing vector normal to the horizon has important consequences 
for the physics of a black hole. In particular, it can be used (with suitable additional assump- 
tions) to prove that the surface gravity of the horizon is a constant [39]. Thermodynamics 
of black holes identifies the surface gravity with temperature, so existence of a Killing vector 
normal to the horizon can be thought of as the condition that the black hole is in equilib- 
rium. Moreover, for stationary black holes with compact horizons, existence of such vectors 
can often be established using Hawking's theorem [5TH I5T~]. This reinforces the intuition that 
stationary metrics with compact event horizons correspond to equilibrium black holes. 



The second-order metric (75) is stationary but the horizon (86) is not compact. On phys- 
ical grounds, we do not expect our black hole (or, rather, black brane) to be in equilibrium; 
indeed, it is dual to an expressly nonequilibrium, entropy-producing process — propagation of 
a shock wave. This means that the condition that the horizon is compact, in the statement 
of Hawking's theorem, is not a mere technicality. Indeed, recall that the time-independence 



of the metric ( 75 ) comes about because we can boost to the frame where the shock wave is 
at rest. This will not be possible, for instance, in global coordinates, where the boundary of 
AdS is a sphere. In that case, the horizon is compact but no longer time-independent: we 
expect it to contain a solitary spherical wave that propagates, say, from north to south and 
eventually equilibrates near the south pole. 



Even though the metric ( 75 ) does not satisfy the compactness condition, and so Hawking's 
theorem cannot be used to establish existence of Killing vector normal to the horizon, it is 
interesting to ask if such a vector nevertheless exists. That would mean that the metric has 
a hidden symmetry that we have up to now failed to detect. We therefore ask if there is a 

X M (r,x) = xto)H*)] + k M (r,x) (94) 



deformation of (93) 
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that satisfies the Killing equation 



V M X N + V N X M 



(95) 



up to and including O(m^) terms. Note that (94) has the same structure as our earlier 



expansions in in particular, it starts with the equilibrium result now considered as a 
function of a slowly-varying parameter. Since X(o) remains a Killing vector to the first order 
and, up to a trivial rescaling, is the only such vector that is null at the horizon, we may 
assume that k M in (94) is 0{u^). This allows us to neglect derivatives of k M with respect 
to x: 



k M (r,x) 



Mi 



(96) 



Substituting (94) into (95) and expanding to the second order in u^, we find that vanishing 
of the (yy) and (zz) components of the Killing equation requires k r = 0, but then the (rr) 
component is explicitly nonzero: 



2u 
r 



2 r 3 T r' T ^0 



(97) 



We conclude that there is no Killing vector of the form (94) at the second order. 



4.4 Expansion of the horizon 



The Killing equation (95) is equivalent to 



£ x 9mn — . 



(98) 



where C x is the Lie derivative along a putative Killing vector. The absence of a solution 
at the second order means that there is no symmetry of the entire spacetime whose Killing 
vector would be normal to the horizon. It is then natural to ask if there is a symmetry of 



the horizon alone that may do instead. This motivates one to consider, instead of (98), the 
quantity 



h 



ab 



£y9 



ab 



(99) 



where g a b is the induced metric on the horizon, and V is a tangent vector constructed as 
follows. Consider the normal (87). Since V M = g MN VN is tangent to the horizon, it can be 



written as the pushforward of some vector V a that lives in the tangent space (and has one 
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fewer component than V 



a M _ 
M = <^C_ 



(100) 



Here, x are coordinates in the spacetime, and £ a are coordinates on the horizon. V is the 



vector used in (99). With the choice £ a = (t,x,y,z), the relation (100) becomes simple: 
V a = V a for all a. 

In general, vanishing of h a b is perhaps the broadest sense in which a horizon can be 
considered equilibrium or non-expanding. As such, this criterion has appeared in a number 
of recent papers (for a review, see [B2j)[^| We will now see that it is not satisfied in our 
present case. 

Although, as we will see, h a b is nonzero already at the second order in Uoo, for use in the 
subsequent discussion, we compute it to the third order. For this, we need g a b and V a only 
to the second order (all the derivatives in the computation of h a t, are taken in the tangent 
space, so each produces an additional power of Woo). The nonzero components of the induced 
metric are 



g tt {x) = r 2 H (x)u 2 (x) + H a (r H ,x) , 

9tx(x) = -r 2 H (x)u°(x)u(x) + H tx (r H , x) + u°r' H (x) , 

9xx{x) = r 2 H (x)ul{x) + H xx (r H , x) - 2ur' H (x) , 

9yy( x ) = 9zz(x) = r 2 H (x) - -s T (r H ,x) 



(101) 
(102) 
(103) 
(104) 



Using the expressions (81)-(83) and working to the second order, we simplify the first three 
entries into 



9tt(x) 

9tx(x) 
9tx(x) 



r 2 H (x)u 2 {x) + ^s T (r H ,x) 



— rff(x)u° (x)u(x) 



-s T (rH,x 



(105) 

(106) 
(107) 



In addition, to this order, 



s T (r H ,x) = (\n2 + 7rT (0) u' {1) (x) 



(108) 



3 We thank A. Ashtekar for bringing [S5] to our attention. 
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and 



V a (x) 



(u°(x),m(x),0,0) 



(109) 



Note that there in no standalone correction for V a : the only second-order terms are those 
contained in the slowly varying parameters u and u°. 



Computing the Lie derivative (99), we obtain the following third-order result: 



h a b(x) = q(x) 



( u 2 (x) -u(x)u°(x) \ 
-u(x)vP(x) [u°(x)) 2 



V 



o 
o 








(110) 



I J 



where 



q{x) 



4 1 

-r 2 H (x)u'(x) + — s ' T (r H ,x) . 



:ni) 



The leading nonvanishing contribution to h a b is of the second order in u^] in (110), we have 
such leading term plus the next-order correction. 
The leading order h a b satisfies the curious relation 



h ab( x ) = 2vr 2 T ( 2 0)( r a6 



;il2) 



where 



with 



a 



fJbV 



-P a P/{d a u p + d p u a ) - -P^d a u a 



P, 



fj,v 



;ii3) 



(114) 



determines the first dissipative correction to the stress tensor of a conformal fluid (in the flat 



4-dimensional spacetime). Thus, (112) relates the expansion of the horizon, as measured by 



the Lie derivative (99), to entropy production in the dual fluid. 



At the next order, the simple relation of h a b to the stress tensor of the dual fluid breaks 
down. It is worth noting, though, that the definition of entropy production at that order is 
already ambiguous [39J. 

We conclude this section by extracting a scalar measure of the horizon's expansion. For 
that, we need to define an inverse g ab of the induced metric Since the vector V, eq. 



(109), is a null eigenvector of g a b, such an inverse is uniquely defined only within the subspace 



orthogonal to V. However, since V is also a null eigenvector of the expansion tensor (110) 
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the scalars we define below will be insensitive to this ambiguity. 

Consider then the projector on the subspace orthogonal to V. To the second order 
inclusive, it has the form 



P%(x) = 



/[«(*)] 



( u 2 (x) -u°(x)u(x) \ 

-u°(x)u(x) [u°(x)} 2 

f[u(x)} 

f[u(x)} J 



where f[u] = (u ) 2 + u 2 . Define the inverse metric by the relation 



(115) 



g ab (x)g bc (x) = P a c (x) • 



(116) 



A convenient diagonal solution is 



g a \x) = 



r H(x)f[u(x)} 



( 1 \ 

10 

f[u(x)] 

V f[u(x)\ J 



(117) 



The requisite scalars are the eigenvalues of the matrix h a c = g h bc . The latter are Xi(x), 
and — |Ai(x) (the last one with multiplicity 2), where 



r%{x) 



(118) 



which measures the expansion of the horizon. 



5 Corrugation (in) stability and sound waves 
5.1 Matching conditions 

Shock waves can be unstable under oscillations of the shock interface. For small perturba- 
tions, stability can be tested by a linearized approximation where the existence of exponen- 
tially growing modes corresponds to an unstable situation. Even if a shock is stable, there 
is still the interesting possibility of a mode which is not damped. Once such a perturbation 
is established, it does not decay in time. This usually describes emission of a wave (e.g., 
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1 

supersonic 



2 

subsonic 




Figure 5: A sketch of the corrugation perturbation. In the supersonic part of the flow there 
is no sound or shear (entropy) wave propagating. The arrow indicates a normal to the 
corrugated shock front. Conservation of the stress energy tensor across the surface is given 



by (119) 



sound) by the shock, the energy for such emission being provided by the shock itself. In this 
section we make the linear analysis for the relativistic shock. As a first step however we need 
to establish the matching conditions for a shock of generic shape. If the shock is described 
by a surface whose normal is everywhere space like and denoted by n M with n 2 = 1 then the 
matching conditions are simply 

n,Tr = n^Tf. (119) 

This is just the continuity of energy and momentum flow through the shock. For the energy 
momentum tensor we take the ideal fluid case 

= T 4 {rf + 4uV) , (120) 

and therefore we get 

T\ {n v + 4(n ■ u x )u\) = T 2 4 (n" + 4(n ■ u 2 )u v 2 ) . (121) 

This equation can be thought as a way to determine T 2 ,«2 from 7\ and u\ (or vice versa). 
In order to do so we project the equation along n M , u\ and u v 2 . This gives rise to three scalar 
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Figure 6: A sketch of the shear or entropy wave. 



equations which can be recast as 



1 + {n- Ul f + {n-u 2 f -8{n- Ul ) 2 (n-u 2 ) 2 = 0, (122) 
ui ■ u 2 + 2{n ■ u\){n ■ u 2 ) = 0, (123) 
if (1 + 4(n • Ul ) 2 ) = Tt (1 + 4(n • u 2 f) . (124) 



The first equation gives 



, 1 + (n ■ ui) 2 . „„. 

( "^ ) = '/ >.^-i (125) 



and the other two can then be used to trivially obtain («i ■ u 2 ) and T2. Notice that we need 
(n-Ux) > This is equivalent to saying that the velocity perpendicular to the interface is 
Vi > |, a constraint that we already knew. Now we would like to obtain the vector u 2 . It is 
clear that u v 2 is a linear combination of u v x and n y and it is not difficult to see that 

u , = 3(n-u 1 )u» + (l-2(n-u 1 ) 2 )n» 
2 v/l + ^-u^n-Hi) 2 -! ' 



5.2 Linearized analysis 

Consider a flat stationary shock wave perpendicular to direction x, namely = (0, 1, 0, 0). 
If a perturbation that deforms the shock wave is produced, it generates waves that move at 
most at the speed of sound with respect to the fluid and therefore do not propagate into the 
supersonic region. So, at the linear level we need to solve 

d^TJf 1 ' = 0, subsonic side, (127) 
5n u (Tf - TD = n^Tf, interface. (128) 
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If we make a ripple on the surface given bjj^Ja; = ( = C,Qt %kyy iut then, to the first order in £, 
the variation in the normal is 

5n» = -i((u,Q,k y ,Q) . (129) 



We also propose that 



with 



kpSTf = . (131) 

There are two types of waves that solve these equations: shear waves and sound waves given 
by 

6T sh = 0, Svt h = - V ^Svf, kl h = -, 

ui y v 2 

6^ K-p CT- (1-*-^-, (132) 

T 2 uj — k x v 2 T 2 ui — k x v 2 

2 

k 2 x + k 2 y - uj 2 = 2 (cj - k x v 2 ) 2 . 
l v 2 

These are shear and sound waves boosted with velocity v 2 . In a fluid at rest a shear wave is 
any time independent perturbation of the form Sve lhx with k ■ x = 0, which in the absence 
of viscosity is a solution of the equations of motion. Notice also that u and k y are the same 
for both waves but the values of k x can be and actually are different. Finally matching at 
the interface we get the equations 

5vl h + 6v% = -hc co(l + 3v 2 2 ) , (133) 

2 



,sh , x „,s _ - u /■ 1 ~ 3t; 2 

3t> 2 

5T S 4 i( uv 2 



K+H = iA *Co-g^, (134) 

(135) 



T 2 3 1 - v\ ' 

Therefore we have six equations for six amplitudes (5v^ ,8v^ ,5v x ,SVy ,ST S ,(o) ■ The charac- 
teristic polynomial has to vanish, which gives a relation between u and k y : 



w = 7!Vr^' (136) 

6 Of course the coordinates are real, we do the standard trick of defining them complex and taking the 
real part at the end. 
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plus another branch where u has the opposite sign compared to k y . Furthermore the disper- 
sion relations determine: 

k x = -- -^uj±uj, (137) 

1 - 3vi 



for the sound wave and 



kf 1 = ~ , (138) 

v 2 



for the shear wave. What is quite interesting is that the sound that is produced has just two 
possible values of the ratio ^ meaning that the direction of the outgoing sound is fixed. It 

Ky 

is also important to notice that u is real meaning that there is no instability of the shock. In 
the ideal fluid approximation the oscillations of the shock are not damped and continuously 
produce sound taking energy from the fluid. It would be interesting to extend this to the 
next order where viscosity is included and see how the oscillations are damped. 



6 Discussion 

In this paper we have continued our study of shock waves in strongly coupled plasmas. 
Macroscopically, shocks are associated with a discontinuity of the temperature and pressure 
along some surface. They typically arise in supersonic flows. The appearance of the discon- 
tinuity signals a breakdown of the hydrodynamic approximation except for the case of weak 
shocks where viscosity is sufficient to resolve the discontinuity. 

Using AdS/CFT we have studied the asymptotic exponentially decaying tails of both 
strong and weak shocks for strongly coupled conformal plasmas in arbitrary dimensions. In 
the rest frame of the shock front the asymptotic behavior is encoded in the penetration 
length as a function of the velocity i = ^y. We have numerically determined this curve 
for different dimensions and confirmed the analytic result that asymptotically, for velocities 
close to the speed of light, the penetration length scales as £ ~ 7~ 2 / d . This means that the 
penetration length is an ultraviolet property of the medium. This is seen quite explicitly on 
the gravity side of the correspondence, where at large 7 the spurious pole that controls the 
scaling approaches the boundary of AdS. It would be interesting to examine other theories 
with gravity duals to find the value of the scaling exponent there. In the context of QCD, 
our result suggests that a perturbative calculation may be enough to obtain the penetration 
length of the sound mode for ultrarelativistic shocks. 

For weak shocks, we have also examined the main properties of the analytic solution 
for the bulk black-hole metric originally found in [35]. We have established that there are 
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no non-trivial Killing vectors in this background and therefore conventional surface gravity 
cannot be defined. Even the less strict condition of a non-expanding horizon does not apply 
to our solution. Indeed, we have found that, to the leading order in the gradient expansion, 
the expansion of the horizon is directly proportional to the entropy production by a shock in 
the dual fluid. Another interesting aspect of shock waves is generation of sound. For ideal 
fluid, the corrugation instability calculation shows that the sound is emitted at a definite 
angle to the direction in which the shock propagates. 

There are many interesting aspects of shocks that we have not addressed here. The ques- 
tion whether a full gravity solution, analytical or numerical, exists has not been completely 
answered. Such a solution would resolve the important question of what happens in the 
region where the linear perturbations diverge. These solutions should also be interesting 
objects in gravity as we expect that they do not have constant surface gravity at the hori- 
zon. Another direction of study would be to understand what happens to the corrugation 
(in) stability when viscosity is taken into account. The inclusion of viscosity introduces a new 
length scale, the resolution length of the shock, and this means that multiple scale analysis 
has to be used. The physical expectation is that the sound emitted by the shock will be 
damped by the viscous terms. 

The recent calculation of thermodynamic functions and transport coefficients for super- 
fluids in AdS/CFT, for example in [53j,[54], (55], also raises the question of existence of shocks 
in superfluids. Shock waves in superfluids have been studied for example in [56]. Two intrigu- 
ing possibilities arise. First, it is possible to have a discontinuity in the superfluid component 
of the velocity of the fluid while the normal component remains continuous. Second, since 
the two sides of the shock have different temperatures it is possible to have different phases 
on the two sides of the shock. We leave these interesting possibilities for future research. 
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A Derivation of the equations of motion 



In this appendix we describe the method to derive the equations of motion for the sound 
channel. The background metric is given by 



7 2 1 r // 7 v Z 2 ( i 1 7 V \ 

as = —ri av dx^dx H — r at , = — dx - = H 



(139) 



and the perturbation by 

iqx 

ds 2 = —j- (h 00 dt 2 + H n dx 2 + 2H 01 dtdx + H(dx 2 d _ 2 )) . (140) 
Einstein's equations 

Gmn = Rmn + dguN = , (141) 
give seven independent equations namely the zz, zt, zx, tt, xx, xt, Xd-i%d-2 components of 



(141). We form the linear combination 

C = A MN G MN , (142) 

and we can choose four components of A, A zz , A zx , A xx , A Xd _ 2Xd _ 2 to eliminate Hq%, H u and 
their derivatives. The remaining two constants (the third is related to an overall normal- 
ization and is not relevant) are not sufficient to eliminate either of the remaining functions. 
However a change to 

H 00 (z) = Z(z) - g(z)H(z) , (143) 



with 



«M = i + ^^£. ( 144 ) 

2 1 — v z 



allows one to eliminate H(z) and write a decoupled equation for Z(z) (28). Then the other 
perturbations are determined from Einstein's equations where Z(z) generically appears as a 
source term for the remaining fluctuations. We do not give here the equations satisfied by 
the other perturbations but an interested reader can find the relevant equations for d = 4 in 
appendix A of 
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B Regularity at the horizon 



The boundary conditions imposed for the perturbations on the the horizon (49) are such 
that the perturbations diverge for positive imaginary values of q. This is often the case 
for perturbations around a black hole. In this appendix we wish to examine whether this 
divergence is an artifact of the coordinate system we use or an effect with a physical origin. 
Let us restrict ourselves without loss of generality to the case of four dimensional plasmas 
and five dimensional black holes that was originally examined in [35]. One can construct 
several covariant quantities that give information about the geometry. Ideally we would like 
to know all the eigenvalues of R a ^s- However this goes beyond the scope of this appendix. 
Out of the 25 invariants that can be built from the eigenvalues of the Riemann tensor in five 
dimensions we examine only one. Since both R^ U R^ V and R are regular from the equations 
of motion we will restrict ourselves to calculating R a j3ryS Rj 3 ' yS which will give us a first insight 
into regularity. Evaluating R a /3 ^ s Rj 3 ' y5 to the first order in the perturbations ( 26 ) and using 



the equations of motion which can be found in appendix (B) of [35] we find that 

r„n» = 8(54 J teS) - 28 8 i— H(r) + ra j 1 ," 

z h z h {z 4 - 3^(1 - v 2 )) 



{ 288(1 _^(3^(3» a -5)) + 4gg2 



z 6 



zl^-Sziil-v*)) 



e- iqx Z(z) 



(145) 



The expression ( 145 ) has been evaluated on shell, meaning that all derivatives higher than the 
second have been substituted using Einstein's equations. We now have to use the behavior 
of the perturbations close to the horizon to examine whether the invariant quantity diverges. 
Using table from [33], which for reference is reproduced here in table (|B]) we find that close 
to the horizon 

Thus, there is a curvature singularity at the horizon. A way out of this would be to suppose 
that as the metric diverges close to the horizon the non linearities become important and 
smooth out the curvature. Such a line of thought can only be pursued numerically and we 
leave it for future work. 

The result of this appendix may appear to be in contradiction to the construction of a 
smooth metric describing weak shock waves in [33] . However this is not true. The reason 
is that, in the Eddington-Finkelstein coordinates used in [33], the location of the potential 
divergence at the horizon is pushed to the past infinity. To see that, let us use the coordinate 
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future singularity 



future horizon 



bifurcation point 




AdS boundary 



past singularity 



Figure 7: A rough sketch of the Penrose diagram for a boosted black hole. Only the rele- 
vant parts are shown. Starting from a patch in the boundary and keeping x, t constant in 
Eddington-Finkelstein coordinates we construct tube shown in cyan. On the other hand if 
we keep x, t constant in Poincare or Schwarzschild coordinates we obtain the tube shown in 
red. The potential divergences arise in the region close to the bifurcation point. In order to 
explore the same problematic region in Eddington-Finkelstein coordinates one has to take a 
spatial and temporal infinity limit . We believe that the resolution of the divergences will 
appear only in the full non-linear solution. 

change from the Poincare coordinates to the Eddington-Finkelstein ones close to the horizon: 



x = x — 



log(z h - z) 



(147) 



This changes the boundary conditions at the horizon from (49) to 



Z ~ const. , Z ~ (zh — z) zVi-v 2 . 



(148) 



That is the infalling condition becomes a constant and the outgoing becomes a non-analytic 
function for generic values of q. This has been used in [35] to argue that the boundary 



value problem set by (28), (49) is well posed. It should be noted here that the location of 
horizon in Poincare coordinates z — z^ with x,t finite translates into z = Zh with x,t going 
to — oo. That is there is no contradiction with the construction in the Eddington-Finkelstein 
coordinates since we are not examining the same geometric region. A schematic of this 
argument can be found in figure ([7]). 

Having calculated the invariant ( |145 ), we can also examine its behavior close to the pole 



of the sound propagation equation (28). From the analytical solution close to the pole, we 



know that it is a regular point. So, we do not expect any divergence, even though some of 



the coefficients in (145) have second or first-order poles at z — Zf. Using the solution around 
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z = z f (35), (38 



0,(39) we find that the curvature invariant is indeed finite at this point. 



Z(r) 



r — > oo 



,-4 



r -> r h 



(r — rh) 4r hV(i-v' 2 ) 



r —7- Tf 

Co + Ci(r - rj) + C 2 (r - r^ 2 + 



#oo(r) 



,-4 



3d 2 -2 



(r-r h ) 4r fe V(i-- 2 ) ) v < yjl 



,(o) 



(r - rft) 4r hV(i-" 2 ) ? v < 



h + hi(r - r } ) + 



2iv{-iqv+2(l-2v 2 )Vl^v 2 -q 2 {l-v 2 ) :i / 2 ) 
(3v 2 -2)(qv-Ai^l-v 2 ) 

qv — 

— 'I 7 — 

•{r — rh) 4r ^ (1 "" 2) , f < 



#(r) 



C 2 \A" - r ft ,w > 



#oi(r) 



m , /,(0D (r _ r/) ^7 



/ip^~(r — r/) + ■ ■ ■ 



#11 (r) 



Table 1: The behavior of the perturbations close to the boundary, the horizon and the pole 



of the equation (28) is given. To connect with the notation used in the present paper, we 
i r h = i. Here, r f 



have r 



j_ 

2 / 



i 



rh 



is the inverse of the location of the 



3I/4 (l_„2)l/4 

pole of (|28l). The boundary conditions for the perturbations at the boundary r — > oo are 



that the metric is unchanged from the Minkowski metric. At the horizon the condition is 
that the asymptotic behavior corresponds to an infalling graviton in the unboosted black 
hole metric. In the linearized approximation the overall scaling factor does not appear and 
we only show the power law behavior. The normalization of the perturbations close to the 
horizon are relative to the normalization of Z(r). The normalization of Z(r) is taken to be 



1 for the factor multiplying (r — r^) 4r hV(i-v 2 ) . p or v > 
not computed. 



the relative coefficients C\ 2 are 
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